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Abstract 

A double-folding method is used to calculate the nuclear and Coulomb interaction between two 
deformed nuclei with arbitrary orientations. A simplified Skryme-type interaction is adopted. 
The contributions of nuclear interaction and Coulomb interaction due to the deformation and 
orientation of the nuclei are evaluated for the driving potential used in the description of heavy-ion 
fusion reaction. So far there is no satisfactory theory to describe the evolution of the dynamical 
nuclear deformation and orientations during the heavy- ion fusion process. Our results estimated 
the magnitude of above effects. 
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I. INTRODUCTION 



The activity of the study on the synthesis of super-heavy elements is still maintained in 
both experimentally and theoretically On the experimental branch, S. Hofmann et. al 
in GSI, Darmstadt, performed experiments on the synthesis and identification of the nuclei 
2 111 and 112 in order to confirm their previous results obtained in the middle of 90s of 
last century (3, [J. Furthermore, several additional decay chains from the reactions 64 iVi+ 208 
Bi ^ 273 111* and 70 Zn + 208 Pb ^ 278 112* were also measured. The joint IUPAC-IUPAP 
Working Party (IWP) has confirmed the discovery of element with atomic number 110, which 
is named as darmstadtium (Ds), recently the new element with atomic number 111 has also 
been proposed by IWP to be named as roentgenium (Rg). Experiments on the synthesis 
of new elements with atomic numbers 115 as well as 113 in the reaction 2A3 Am + 48 Ca 
were carried out at the U400 cyclotron in Dubna recently they also reported the results 
of excitation-function measurements for the 244 Pu + 48 Ca fusion-evaporation reactions for 
element 114 and the synthesis of new isotopes of element 116 with the 245 CVn+ 48 Ca reaction 

Q- 

On the theoretical branch, the physics on the more complicated dynamical process to 
super-heavy elements has been paid more attention [(J and investigated by several groups 
under different mechanisms, for example, the di-nuclear concept (see the^ recent works in 
[AM and the references therein), the fluctuation-dissipation model 10, 

nn j nn 

121,1121], as well as the macroscopic dynamical model |14L 115]. 



nucleon collectivization 



the concept of 



In the di-nuclear system (DNS) concept [llSIS) 16, l?l 3> the fusion process is considered 



as the evolution of a di-nuclear system caused by the transfer of nucleons from the light 
nucleus to the heavy one. The nucleon transfer process is described in Ref. [8( by solving 
the master equation numerically. It is found that the fusion probability of the compound 
nucleus is very sensitive to the specific form of the driving potential. In Ref. Q], the Coulomb 
interaction potential of deformed nuclei with a tip-tip orientation is considered. However, 
spherical nuclei were adopted in calculating the nuclear interactions since it is thought that 
the nuclear interaction does not depend so strongly on the deformation of the nuclei as the 
Coulomb interaction due to the short range characteristics of the nuclear force. Although 
some reasonable results, such as the optimal excitation energies, the residual cross sections 
of super-heavy compound nuclei, were obtained for different heavy-ion fusion reactions, the 
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reliability has to be checked. 

Presently a double-folding method is developed to calculate the nuclear and Coulomb 
interactions between the two deformed nuclei with arbitrary orientations. Here we consider 
the ground state deformations of the nuclei for all possible combinations of the DNS of a 
certain reaction. In principle, the deformed nuclei can have different relative orientations 
which supply quite different conditions for fusion. Some averaging over the orientations of 
the nuclei has to be carried out at least in the entrance channel. The deformation and the 
orientation evolutions are difficult to be described, which have not yet been investigated 
very well by any model so far. Nevertheless, it is important to bear in mind : what are 
the magnitudes that the deformation of nuclei contributes to the nuclear and Coulomb 
interactions, respectively, and to explore how and to which extent the orientations contribute. 
These investigations will give a direction for further improvement. 

The paper is arranged as follows. In the next section, the treatment of the nuclear and 
Coulomb potentials is introduced. We present the calculated results and the corresponding 
discussions in Section III, where the interaction potentials between different deformed nuclei 
and their dependence on orientations as well as the driving potentials used in the DNS model 
for different fragmentations in reactions leading to 272 Ds. Finally, Section IV gives a brief 
conclusion and outlook. 



II. TREATMENT OF DRIVING POTENTIALS FOR ORIENTATED DE- 
FORMED NUCLEI OF DNS 

For a dinuclear system, the local excitation energy is defined as follows, 

^Urel int 

where E* is the intrinsic excitation energy of the dinuclear system converted from the relative 
kinetic energy loss, M is the corresponding intrinsic spin due to the dissipation of relative 
angular momentum J. J re \ and Ji n t are the relative and intrinsic moments of inertia re- 
spectively. XJ(A\,A-2) is the driving potential energy responsible for the nucleon transfer in 
the DNS model, and is written down as, 



3 



U(Ai, A 2 , R) = U ld+ sc{Ai)+U ld+ sc(A 2 )-U ld+ so(Acn)+Uc(A 1 , A 2 , R)+U n (A 1 , A 2 , R), 

(2) 



where Ai, A 2 , and Aqn represent the mass numbers of the two nuclei and the corresponding 
compound nucleus, respectively, we have A\ + A 2 = Aqn- In the DNS model, the driving 
potential is normally given as a function of rj = (A\ — A 2 )/Acn- The first three parts of the 
right hand side of the equation are calculated from the Liquid-Drop model plus the shell and 
pairing corrections Uc{A\, A 2 , R) and Un(Ai, A 2 , R) are the corresponding Coulomb 

and nuclear potential energies between the nuclei and depend on the fragmentation of the 
dinuclear system, on the internuclear distance R and on the orientation and deformation of 
the nuclei. They could be calculated by different methods. In the present work, we calculate 
them by using the double-folding method 



where px( r x) and p 2 (r 2 ) are the density distribution of 1 and 2 nucleons in a dinuclear system, 
v(rx — r 2 ) is the corresponding interaction between the two points. For the nuclear part Un 
we use densities with a smooth falling off at the surface (see later) and constant densities for 
the Coulomb interaction. The long-range Coulomb interaction is not sensitive to the density 
at the surface which allows to simplify the calculations. Therefore, we write the Coulomb 
interaction as follows, 



where R is the vector between the two centers of the nuclei ("T" and "P") as illustrated 
in Fig^ The charge densities are set as p\ = |^ and p\ = f^f where Z 1>2 and fi 12 are 
the proton numbers and the volumes of the two nuclei, respectively. The symmetry axes 
(ax and a 2 ) of the two deformed nuclei and the z-axis are assumed to be in the same plane. 
7i and 72 are the corresponding angles between the symmetric axes and i-axis, i.e., which 
represent the different orientations of the two nuclei, while ot\ and a 2 are the angles between 
arbitrary vectors r 12 and the symmetric axes Si and a 2 , respectively. The distance between 
the two points is given by 




(3) 




(4) 




(5) 
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FIG. 1: Schematic presentation of the orientation of two axially quadrupolly deformed nuclei. 
It is easy to find the following relations, 

(ri - r 2 ) 2 = r 2 + r 2 - 2rir 2 (sin6'i sin# 2 cos(0i - <f> 2 ) + cos#i cos# 2 ), (6) 

(ri - r 2 ) • R = (n cos 9 1 - r 2 cos 9 2 )R. (7) 

where 9i >2 and 0i j2 are the angles of r ii2 with respect to the coordinates (x, y, z) and (xt, 
yf, zf), respectively. 

The upper and lower limits of r 12 , 12 , and 4>i,2 are 

ri,2- (0,% )2 )); 0i l2 :(O,7r); lj2 : (0, 2vr), (8) 

where $ft(aii) and 3?(o; 2 ) describe the nuclear surface with quadrupole deformations. 

to(a i ) = R Qi (l + ? 2 Y 20 (a i )). (9) 



Here Roi are the spherical radii of the two nuclei which preserve their fixed volumes. Y 2 o(a) = 
(5/47r) 1/,2 P2(cos a) = (5/47r) 1//2 (3 cos 2 a — 1)/2 is spherical harmonics and the axial symmetry 
is preserved. The /3| is the quadrupole deformation parameter of z-nucleus taken from 
Ref . . It is easy to write down the expressions for aii and a 2 as 

cos«i = a 1 • = sin^x cos0i sin 7! + cos^ cos 71, (10) 

and 

cos a 2 = a 2 ■ ^(0:2) = sin 9 2 cos <p 2 sin 72 + cos 8 2 cos 72. (11) 



For the nuclear potential, following the work by Adamian et al. [18], we adopt the 
Skyrme-type interaction without considering the momentum and spin dependence, in which 
a zero-range treatment of the effective interaction S(ri—r 2 ) is assumed. The nuclear potential 
is obtained in the sudden approximation jl8^ . 



U N (R) = C { m ex (l ' pl(r)p 2 (r-K)dr 
Poo J 

+ J pi(r)p 2 (r - R)dr) + J pi(r)p 2 (r - R)dr} (12) 



with 



f. -f +r ^-z^-z* (1 * 

r m,ex — Jin,ex ~r J i n ex A A ■ \ 10 J 

Here Ni j2 and Zi j2 are the neutron and proton numbers of the two nuclei respectively. 
Obviously the isospin effect of the nucleon-nucleon interaction is considered here though the 
relative influence is small. The parameters Cq = 300 MeV-fm 3 , /j n = 0.09, f ex = —2.59, 
//„ = 0.42, f' ex = 0.54, and p o = 0.17fm~ 3 are also used in this work. The functions pi 
and P2 are two-parameter Woods-Saxon density distributions (now we set the center of the 
"P" -nucleus at the coordinate origin and ri = r) 

91 W = l + exp((r-^ 1 (a 1 ))/a pi ) ( U) 

and 

P2W = l + exp((|r-R|-K 2 («2))/a P2 )' (15) 
The parameters a pi and a P2 represent the diffuseness of the two nuclei, respectively. Whereas 



cos«i is given in Eq. (jlUJl . we use the following formula with |r — R| = \Jr 2 + R 2 — 2rR cos 9 
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r-R| v ; 

r (sin 9 cos <j> sin 72 + cos 9 cos 72) — R cos 72 
r 2 + R 2 -2rRcos9 ' 

We directly calculate the six- and three-dimensional integrals in Eqs. (J1J and (fE?j) nu- 
merically. For Eq. ()12|). a truncation parameter r cut for the upper limit of r is introduced 
due to the long tails of the nuclear densities expressed in Eqs. ()14|) and (fTo^) . For each 
mass asymmetry we calculated the sum of the Coulomb and nuclear potential energies as a 
function of the internuclear distance R and took the potential at the minimum in R which 
is shorter than Rqb (Coulomb-barrier saddle point) as the driving potential of the DNS 
model. 



III. NUMERICAL RESULTS 



In this paper, the nuclear and Coulomb interaction for the DNS of the reaction 64 Ni + 208 
Pb —>- 272 Ds is studied by taking the nuclear deformations and the corresponding orientations 
into account. For simplicity, the diffuseness parameters a Pl and chosen 
0.6fm, which is a little bit larger than those in Ref. Q|. Furthermore, roi = = 1.2fm is 
used. The parameter r cut = 25fm for the radial integration of the nuclear potential of the 
deformed nucleus in Eq. (I12|) . is taken for an adequate precision. 

Fig|21 (a) and (b) show the nuclear interaction potentials of two sets of projectile-target 
combinations, namely 2S Na + 244 Es and 74 Zn + 198 H g, to form the same compound nucleus 
272 Ds as a function of distance R between the centers of the two nuclei. The corresponding 
nucleus-nucleus potentials including both the nuclear and Coulomb interactions are given 
in Fig|21 (c) and (d). In Fig|21 (a) and (c), both nuclei are prolately deformed, 28 Na with 
p 2 = 0.257 and 244 Es with f3 2 = 0.234, respectively, while in FigO(b) and (d), u Zn is prolate 
and 198 Hg oblate with /3 2 — 0.125 and —0.112, respectively. The system 74 Zn+ 198 Hg is more 
mass-symmetric, i.e., it has a smaller \rj\ than the system 28 iVa + 2AA Es, and thus a higher 
Coulomb potential energy. In each panel, different orientations for the two systems, i.e., tip- 
tip, tip-belly and belly-belly orientations are investigated, an illustration is shown in (c) plot. 
When both (3 2 values are positive in (a) and (c), the angles (71,72) = (0°, 180°), (0°,90°), 
and (90°, 90°) are the corresponding ones for the tip-tip, tip-belly, and belly-belly cases, 



respectively, while for the case of f3\ > and /3| < in cases (b) and (d), (71, 72) = (0°, 90°), 
(0°, 0°), and (90°, 0°), respectively. The two nuclei become more compact with a belly-belly 
orientation in contrast to the tip-tip one, and the minimum of the potential energy for a 
belly-belly orientation is at a smaller R than that of the tip-tip case. When the orientation 
changes from the tip-tip type to the belly-belly one, the minima of the nuclear potentials in 
(a) and (b) behave differently as those of the total potentials shown in (c) and (d), i.e., the 
minima of the nuclear interaction go down while the minima of the total interaction increase. 
The reason for the decrease from tip-belly to belly-belly in (c) is that the increase of the 
Coulomb interaction energy is smaller than the decrease of the nuclear interaction energy. 
Defining a distance between the surface of the two nuclei, for example, for the tip-tip case, 
AR = R min - ($r 9 + $2° ng )i while for the belly-belly case, AR = R min - (Uf ort + ^f ort ) ( 
^iong,short re p resen t ^he long and short axes of the deformed nucleus i, respectively), we find 
that AR changes a little for different orientations. When \r]\ decreases from 1 to 0, the value 
of AR increases due to a larger repulsive Coulomb force, which can be seen more clearly in 
FigE3 Therefore, the effect of the mass asymmetry and the orientation of the DNS on the 
driving potential can be analyzed from these results. 

FiglHl shows the potentials at the minimum of Un + Uq illustrated in Fig|2]for the above 
two combinations as a function of the orientation, the orientation is chosen in a way that 
keeps 7l + 72 = 180° for the system 28 Na + 244 Es and 72 + 71 = 90° for 7A Zn + 198 Hg. On 
the left hand side, 71 goes from 180° to 0° and 72 from 0° to 180°, on the right hand side, 
71 is chosen from 0° to —180° and 72 from 90° to 270° in order to obtain similar trends of 
the variation of potentials as a function of the orientation of the two nuclei as on the left 
hand side, in the two cases, both of the orientation changes from the tip-tip orientation to 
the belly-belly one and finally back to the tip-tip orientation (the orientation of the nuclei 
is shown in the lower-left plot of FigOJ). With the changing of the orientations of the two 
nuclei, the nuclear potentials (upper panels) form a valley while the Coulomb potentials 
(middle panels) attain a peak value for the tip-tip orientation. The summation of the two 
contributions shown in the bottom panels is similar in shape to the Coulomb potential but 
the change with angle is gentler. 

Fig0] displays the driving potentials in Eq. (J2J) for different orientations. In the upper 
panel of the figure we fixed 71 and 72 to 0° or 90°, while in the lower panel, the results for 
the tip-tip and belly-belly orientations are shown. The curves in FigHJwere calculated by 
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FIG. 2: The nuclear (in (a) and (b)) and the nuclear+Coulomb potentials (in (c) and (d)) for 
two sets of projectile-target combinations for the same compound nucleus 272 Ds are shown as a 
function of R for different orientations of the two nuclei. 

starting with the initial fragmentation 6A Ni + 208 Pb transferring nucleons in steps of 
one proton or one neutron by searching for the minimum of potential energy. Therefore, 
the potentials are only approximately symmetric with respect to rj = for the tip-tip and 
belly-belly cases while for the cases with orientations of (0°, 0°) and (0°, 90°) in the upper 
panel, this symmetry is lost obviously. From FigEl we find that the driving potential is 
quite sensitive to the choice of orientations of the two nuclei. The driving potential for the 
tip-tip configuration is smaller than that for the belly-belly configuration in the whole range 



of r). This result is in disagreement with that obtained in Ref. |21|, this discrepancy might 
be associated to the different consideration of the fusion process of heavy-ions. 

To evaluate the difference between the different orientations, we show the differences 
between the potential energies of the various cases with respect to the tip-tip case in the 
upper half of FiglHJ where jj beU y- hell y _ jjtip-tip j g snown w jth a line, while the other two cases 
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FIG. 3: The potentials with 7l + 72 = 180° for 28 Na + 244 Es (left panel) and 7l + 72 = 90° for 
Zn+ 198 Hg (right panel). See text for details. 



are shown with different scattered symbols. The differences are peaked in two regions, one 
in \r)\ < 0.5 and the other in \rj\ > 0.5. In each region there exists a large deformation of the 
nuclei, especially when \r]\ is 0.1 ~ 0.4. However, the detailed deformation of the two nuclei in 
each part is different, that is, when > 0.5, the smaller nucleus is almost spherical while the 
larger counterpart is prolately deformed. When \rj\ < 0.5, prolate and oblate deformations 
of the two nuclei occur, for example, for 77 = —0.243, the corresponding configuration is 
103 Mo + 169 Er with a couple of large prolate deformation (fy, /3%) = (0.358,0.304). When 
f] = —0.169, the corresponding 113 Pd+ l59 Gd consists of a negatively (—0.25) deformed U3 Pd 
and a positively (0.28) deformed 159 Gd. The separation distance AR between the surface 
of the two nuclei of the DNS is shown in the lower graph of FigEl Because of the relatively 
large Coulomb potential, AR is stretched when the masses of the two participating nuclei 
become more equal, which has also been shown in FigEl 
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FIG. 4: The driving potentials with different 
orientations for the system 272 Ds. The calcu- 
lation was started with the 64 iVi + 208 Pb frag- 
mentation(see text). 



FIG. 5: (a): the difference between the driving 
potentials of tip-tip case and the other cases, 
(b) : AR as a function of r\ with tip-tip orienta- 
tion. 



For the di-nuclear system 64 Ni + 208 Pb — > 272 Ds, Fig|B] shows the comparison between 
the present driving potential shown by dots and that calculated in Ref. [s| shown by a fine 
line for the tip-tip orientation. In the present calculations, the ground state deformation 
has been taken into account for both the nuclear and Coulomb interactions, and in Ref. ^] 
a parameterized Morse formula |9| without considering the deformation of the nuclei has 
been adopted for the nuclear potential. We find that the two potentials are basically very 
close each other, however, some obvious deviations appear in the relatively large deformed 
regions, for example, around \r]\ ~ 0.2 and \r]\ ~ 0.8. It should be pointed out that a 
deviation also occurs at \rj\ ~ 0. After checking the detailed path of evolution, we find that 
the configurations of the DNS in the two cases are different at this point. For the case with a 
nuclear interaction of spherical nuclei, the combination ( 136 La + 136 /) is preferred, while for 
the one with that of deformed nuclei, a more charge-symmetric combination ( 136 5a+ 136 Xe) 
is taken into account. Obviously, the effect of a large deformation in the deformed region 
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FIG. 6: The comparison of the driving potential by using the ground state deformation and the 
tip-tip orientation for both the nuclear and Coulomb interactions with previous calculations which 
did not taken into account the deformation in the nuclear part of interaction. 

\rj\ ~ 0.2 changes the final path of the evolution near rj = 0. 

IV. CONCLUSION AND OUTLOOK 

A double-folding method used to calculate the nucleus-nucleus potential between de- 
formed nuclei is further developed to improve the driving potential of nuclear fusion in the 
DNS model. By taken into account the nuclear deformation in the nuclear interaction to- 
gether with the Coulomb interaction, the formalism for calculating the driving potential of 
heavy-ion fusion becomes more reasonable. The deformations and orientations of the in- 
teracting nuclei contributing to the nuclear and Coulomb interactions are investigated for 
every fragmentation of the DNS considered. It is natural that the tip-tip orientation has the 
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lowest interaction energy, and may be preferred during the nucleon exchange process. The 
minimum energies of the nucleus-nucleus interaction along the distance between the centers 
of the two nuclei appear at larger distances when mass- asymmetry \rj\ changes from unit to 
zero, which is due to the larger Coulomb force, and is in favor for the quasi-fission process. 
So far a dynamical evolution of the deformation and orientation during the heavy-ion fu- 
sion process is not reasonably treated by the present models to our knowledge, our results 
have estimated the effects of the deformation and orientation of the nuclei, on the driving 
potential. Hopefully it will give a direction for the further investigation and improvement. 

In the next step, we will calculate the fusion probability of various projectile-target com- 
binations with deformed nuclei. Furthermore, when the distance between the surface of 
the two nuclei is elongated the effect of quasi-fission is expected more pronounced. We will 
further consider a two-dimensional potential as functions of the mass asymmetry 77 and the 
internuclear distance R in order to investigate the effect of quasi-fission in the subsequent 
work. 
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